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ABSTRACT 

We searched for 21-cni absorption associated with the Zabs= 1.3647 absorption sys- 
tem toward the gigahertz peaked-spectrum somxe PKS 0237—233 using the Giant 
Metrewave Radio Telescope. A high quahty UVES spectrum shows that C i and C i* 
are detected at this redshift together with C 11*, Mg i, Mg 11, Si 11, Al 11, Fe 11 and 
Mn II. The complex profiles, spread over ~ 300 km s~^, are fitted with 21 Voigt pro- 
file components. None of these components are detected in 21-cm absorption down 
to a detection limit of r(3cr) 3 x lO'^ (or iV(H i)/Ts ^ lO" cm'^ K"!). We de- 
rive log [iV(H i) cm~^]^19.30±0.30 using the Lya absorption line detected in the lUE 
spectrum of the quasar. Mg 11, Si 11 and Al 11 column densities are consistent with near 
solar metallicity and we measure [0/H] > —0.33. Using photoionization models con- 
strained by the fine-structure excitations of C i and C li, and the 21-cm optical depth, 
we show that the C i absorption arises predominantly cither in a warm ionized medium 
(WIM) or warm neutral medium (WNM) in ionization and thermal equilibrium with 
the meta-galactic UV background dominated by QSOs and star forming galaxies. The 
estimated thermal pressure of the gas is of the same order of magnitude over different 
velocity ranges through the absorption profile (2.6 ^ log[P/fc(cm~'^K)] ^ 4.0). The 
gas-phase metallicity corrected for ionization is Z ^ 0.5 Zq with a signature of iron 
co-production elements being under abundant compared to a— process elements by 
~0.5 dex. At ZabsJsl-9, C i absorption is usually associated with H2 absorption arising 
from cold gas in damped Lyman-a systems. This system and the Zabs = 2.139 toward 
Tol 1037—270 are the only two systems known which show that C i absorption can 
also be detected in warm gas provided the metallicity is high enough. Interestingly, 
both the systems are part of unusual concentrations of absorption lines. 
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1 INTRODUCTION 

The diffuse gas in the interstellar medium (ISM) is known 
to exhibit a wide range of physical conditions (temperature, 
density, radiation field and chemical composition). These 
physical quantities are influenced by the local star forma- 
tion, the cosmic ray energy density, photoelectric heating 

Based on observations carried out at the European Southern 
Observatory with UVES in the course of the Large Program " The 
Cosmic Evolution of the IGM" 166.A-0106 on the 8.2 m VLT 
telescopes Kuyen operated at Paranal Observatory; Chile. 



by dust as well as mechanical energy input from both im- 
pulsive disturbances such as supernova explosions (McKee 
& Ostriker, 1977) as well as more steady injection of energy 
in the form of winds from stars (Abbott, 1982). Understand- 
ing these physical conditions and the processes that main- 
tain these conditions are important for our understanding 
of galaxies and their evolution. 

The relative populations of the ground and excited lev- 
els of C I have been used to discuss the particle density, 
the ambient UV radiation field and the temperature of the 
cosmic-microwave background radiation (see Bahcall, Joss 
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& Lynds 1973; Meyer et al., 1986; Songaila ct al. 1994; Go, 
Bechtold & Black, 1997; Roth & Bauer, 1999; Srianand et 
al. 2000; Silva & Viegas, 2002). In the Galactic ISM, fine- 
structure excitation of C I has been used to study the dis- 
tribution of thermal pressure (see Jenkins & Tripp 2001). 
At high redshifts C l absorption is very rare and detected 
in only a handful of systems (~ 15% of the damped Lyman 
a systems (DLAs) at «abs^ 1-9)- Usually, C I absorption 
in these systems is closely associated with low temperature, 
high density gas that often shows associated H2 absorption 
(Ge & Bechtold, 1999; Petitjean et al. 2000; Ledoux, Pe- 
titjean & Srianand 2003; Srianand et al. 2005). The only 
exception is the Zabs= 2.139 system towards Tol 1037—270 
(Srianand & Potitjcau, 2001) that originates from a solar 
metallicity sub-DLA (i.e. systems that show damping wings 
with log N(H l)< 20.3) very much similar to the one we 
shall discuss here in which the absence of detectable H2 ab- 
sorption lines suggests that the gas temperature is probably 
high. 

Like C I, C 11* absorption has been detected iu every 
system where H2 is present (Srianand et al. 2005). However, 
C II* has also been seen in a considerable fraction of DLAs 
(~ 50%) without H2 (Wolfe, Prochaska & Gawiser 2003; 
Wolfe, Gawiser & Prochaska 2003; Srianand et al. 2005). It 
has been found that C 11* absorption is always detected when 
log N(H i) ^ 21 and also in systems with high metalUcity 
irrespective of N{}i i). As C 11* is an important coolant of 
the ISM, the corresponding C 11* column density can be used 
to discuss the cooling rate in the absorbing gas (Wolfe et al. 

2004) . 

The nature and physical conditions in the H I gas can 
also be probed by 21-cm absorption lines. The searches for 
21-cm absorption in DLAs at Zabs^ 2 have mostly resulted 
in null detections with typical spin temperatures > 10^ K 
(see Table 3 of Kanekar Sz Chengalur 2003 and Table 1 of 
Curran et al. 2005). Based on the lack of 21-cm and H2 
absorption in most of the systems it has been inferred that 
most of DLAs actually trace a warm neutral phase (see also 
Petitjean, Srianand & Ledoux 2000). However one of the 
possible uncertainties to interpret the 21 cm data could be 
the unknown covering factor of the H I gas in front of the 
radio source (see Briggs & Wolfe 1983 for a discussion). 

Despite obvious advantages of studying all the above in- 
dicators simultaneously, such studies are possible only for a 
few systems [at 2abs= 2.8110 towards PKS 0528-250 (Car- 
ilh et al., 1996; Srianand & Petitjean 1998), Zabs= 1.944 
towards Q1157-F014 (Wolfe & Briggs 1981; Srianand et al. 

2005) and 2abs= 2.03954 towards PKS 0458-020 (Wolfe et 
al. 1985; Briggs et al 1989; Hcinmiillcr et al. 2005)]. This 
is mainly because only a small fraction of the QSOs with 
interesting DLAs are radio-loud and, the redshifted 21-cm 
frequency range in many cases is not covered by the present- 
day radio telescopes or is affected by strong radio frequency 
interferences. 

The UVES spectrum of the compact source 
PKS 0237-233, that was obtained as part of the ESO-VLT 
large programme "The Cosmic Evolution of the IGM" 
(Bergeron et al. 2004), shows C I, C I* and C II* absorption 
lines spread over ~300 km s~^ in a strong Mg 11 sub-DLA 
system at «abs= 1.365. The wavelength range into which 
the corresponding 21-cm absorption is redshifted can be 
covered by the Giant Metrewave Radio Telescope (GMRT) 



610 MHz band. Here, we present the analysis of this system 
using the high-resolution optical spectrum obtained with 
UVES/VLT, a low-resolution UV spectra obtained with 
lUE and the GMRT radio spectrum covering the redshifted 
21-cm absorption. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 Optical UVES observations 

The optical spectrum of PKS 0237—233 used in this study 
was obtained using the Ultraviolet and Visible Echelle Spec- 
trograph (UVES; Dekker et al. 2000) mounted on the ESQ 
Kueyen 8.2 m telescope at the Paranal observatory in the 
course of the ESO-VLT large programme "The Cosmic Evo- 
lution of the IGM". PKS 0237-233 was observed through 
a 1" slit (with a seeing typically <0.8 ") for ~12 hr with 
central wavelengths adjusted to 3460A and 5800A in the 
blue and red arms, respectively, using dichroic no. 1 and 
for another ~14 hr with central wavelengths at 4370 A and 
86OOA in the blue and red arms with dichroic no. 2. The 
raw data were reduced using the latest version of the UVES 
pipeline (Ballester et al. 2000) which is available as a ded- 
icated context of the MIDAS data reduction software. The 
main function of the pipeline is to perform a precise inter- 
order background subtraction for science frames and master 
flat-fields, and to allow an optimal extraction of the object 
signal, rejecting cosmic ray impacts and performing sky- 
subtraction at the same time. The reduction is checked step 
by step. Wavelengths are corrected to vacuum-heliocentric 
values and individual ID spectra are combined. Air-vacuum 
conversions and heliocentric corrections were done using 
standard conversion equations (Edlcn 1966; Stumpff 1980). 
Addition of individual exposures is performed using a slid- 
ing window and weighting the signal by the errors in each 
pixel. Great care was taken in computing the error spectrum 
while combining the individual exposures. Our final error is 
the quadratic sum of appropriately interpolated weighted 
mean errors and possible errors due to scatter in the in- 
dividual flux measurements. The final combined spectrum 
covers the wavelength range of 3000-10,000A. A typical S/N 
~ 60-80 per pixel (of 0.035 A) is achieved in the whole wave- 
length range of interest. The detailed quantitative descrip- 
tion of data calibration is presented in Aracil et al. (2004) 
and Chand et al. (2004). 

2.2 Radio Observations 

The radio spectrum was obtained on 2004, Feb 9 using 
GMRT. The local oscillator chain was tuned to centre the 
baseband at the corresponding redshifted 21-cm frequency. 
The FX correlator system at GMRT splits the baseband into 
128 channels, yielding a velocity resolution of 3.9 km s~^ for 
our chosen baseband bandwidth of 1 MHz. We observed 
standard flux calibrators 3C48 or 3C147 every two hours 
to correct for the amplitudes and bandpass variations. Since 
the background radio source, PKS 0237—233 is unresolved 
for the GMRT array at the observing frequency, observing 
a separate phase calibrator was not required. A total 8 hrs 
of on-source data were acquired in both the circular polar- 
ization channels RR and LL. 
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Figure 1. Velocity plots of the absorption line system centered at z^bs = 1-3647. The points with error bars are the observations and 
the continuous curve shows the Voigt profile fit to the absorption profiles. Tick marks in each panel indicate the location of individual 
components. We identify 21 distinct components based on the fits to Mg II and Fe II profiles. The distinct velocity ranges discussed in 
detail in the text are indicated in the bottom left panel. 



The radio data were reduced in the standard way using 
the Astronomical Image Processing System (NRAO AIPS 
package). After the initial flagging and calibration, source 
and calibrator data were examined to flag and exclude the 
baselines and timestamps affected by Radio Frequency In- 
terference (RFI). It was noticed that data on 3C 147 was 
not particularly good as compared to 3C 48. Of the 4 time- 
scans (exposures) on 3C 48 and 3C 147, 3C147 contributed 
only to the last one. So we used only 3C 48 for the band- 
pass calibration. Applying these complex gains and band- 



pass, a continuum map of the source was made using line 
free channels. Using this map as a model, self-calibration 
complex gains were determined and applied to all the fre- 
quency channels. The same continuum map was then used to 
subtract the continuum emission (a flux density of 5.1 Jy). 
The continuum subtracted data set was then imaged sepa- 
rately in stokes RR and LL to get three-dimensional (with 
third axis as frequency) data cubes. Spectra at the quasar 
position were extracted from these cubes and compared for 
consistency. The two polarization channels were then com- 
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Figure 2. Excited fine-structure lines of C I and C II plotted on a 
velocity scale with respect to ^aba= 1-3647. Points with errorbars 
are the observed data and the best fit is shown as a continu- 
ous curve. The tick marks indicate the location of distinct Voigt 
profile components. 



bined to derive the final stokes I spectrum, which was then 
shifted to heliocentric frame. 



2.3 lUE spectrum 

An lUE spectrum of PKS 0237-233 was obtained by 
Richard F. Green as part of an observing programme to 
study high-redshift quasars (lUE Program ID: QSDRG and 
lUE image number: LWR10935). The spectrum was ob- 
tained with the large aperture (approximately 10 x 20 arc 
sec) and an exposure time of 25,800 s. The typical spectral 
resolution is ~ 6 A and the signal to noise per pixel ~ 5. We 
used the pipeline calibrated lUE spectra obtained directly 
from Multimission Archive at Space Telescope (MAST). 



3 ^ahs 
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3.1 Neutral and singly ionized species 

The velocity profiles of some of the neutral and singly ionized 
species seen at Zabs= 1.3647 together with a multicomponent 
Voigt profile fit are shown in Fig [1] They are spread over 
~ 300 km s~^. The absorption lines of Mg I, Si ll, Al ll, 
Al III and Fe ll fall in the spectral range with the best signal- 
to-noise ratio. We also detect O lA1302 in the blue part of 
the UVES spectrum that has a lower SNR. Weak Mn ll 
lines are clearly detected in the strongest Mg ll components 
corresponding to velocity ranges B and C (see left bottom 
panel in Fig. 1). We also tentatively detect Cr ll in the same 
components (see Fig. [TJ . 

We detect C l, C l* and C ll* absorption in most of the 
components (see Fig. [2} . Simultaneous Voigt profile fitting 
of all the metal transitions requires 21 independent compo- 
nents. The results are summarised in Table. [T] To facilitate 
the discussion, we divide the velocity profiles into 5 distinct 



Figure 3. Top panel shows the low resolution (~ 5.5 A) lUE spec- 
trum of PKS 0237—233, the fitted continuum (smooth continuous 
curve) and the error spectrum (dashed line). The expected posi- 
tions of the Lyman series lines at 2abs= 1.3647 are marked by 
vertical lines. Single component Voigt profile fits for log N (H l) 
= 19, 19.5 and 20.0 are shown by, respectively, continuous, dot- 
ted and dashed profiles. These fits suggest that the system is a 
sub-DLA with log N{U i) ^ 20.0. 



velocity regions marked as A, B, C, D and E in the bot- 
tom left panel of Fig. [1] The integrated column densities of 
different species within these velocity ranges are reported 
in Table [T] Ni iiA1709, Co iiA1466, Na iAA3303.3,3303.9, 
Ca llAA3934,3969 are not detected, nor are the fine-structure 
lines of O i and Si ii. 



3.2 Constraints on the total Af(H l) 

We used the lUE LW spectrum from the lUE archive to 
derive the total H I column density, A''(H l). We detect a line 
at the expected position of the redshifted Lya absorption 
with rest equivalent width 3.0± l.oA (Fig.|3]). If we assume 
the absorption to be in the damped portion of the curve 
of growth we derive log A'^(H l) — 19.30±0.30. This means 
that N{¥^ l) in individual components cannot be higher than 
this value. If the Lya line is the blend of saturated Lya 
absorptions from individual metal-line components, then the 
total 7V(II l) could be less than this value. 

The derived constraint on A'^(H l) is consistent with 
ASCA & ROSAT observations that do not detect excess 
soft X-ray absorption on top of the galactic absorption (Ya- 
masaki et al. 1998). Thus, the allowed range in A'^(II l) is very 
much consistent with the system being either a sub-DLA or 
a high column density Lyman Limit system. Upper limits on 
N{}i l) in the individual components can be obtained from 
our GMRT spectrum for an assumed spin-temperature (Ts). 
The GMRT spectrum is shown in the upper panel of Fig. ID 
Despite reaching a la limit of 10~^ for r(21 cm) we do not 
detect any absorption at the expected positions of the dif- 
ferent individual components detected in the UV spectrum 
of neutral and/or singly ionized species. 

The neutral hydrogen column density in a velocity in- 
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Table 1. Voigt profile fitting results for individual components: 



b(km/s) 



iV(Mg I) 



JV(Mg II) 



N(0 I) 



iV(Si II) 



JV{Fo II) 



]V(C II*) 



N(C I) 



]V(C I*) 



JV(H I)^ 



Velocity range A 

11.21(0.01) 13.97(0.01) 13.58(0.08) 13.58(0.02) 12.89(0.01) b 11.80(0.19) ^ 12.00 silS.36 

11.38(0.01) 13.47(0.01) 13.62(0.06) 13.56(0.01) 13.04(0.01) b 12.12(0.10) ^ 12.00 :gl8.66 

11.50(0.01) 13.73(0.02) 13.93(0.09) 13.77(0.02) 13.21(0.01) b 12.07(0.12) 12.00 SS18.43 

10.00(0.10) 12.56(0.01) 12.85(0.31) 13.04(0.02) 12.21(0.06) b ^ 11.70 :g 12.00 :gl8.30 

11.86(0.01) 14.26(0.01) 14.23(0.05) 14.16(0.01) 13.56(0.01) b 12.49(0.07) :g 12.30 :gl9.21 



1.364232 2.6(0.2) 
1.364334 5.1(0.2) 
1.364443 3.0(0.2) 
1.364490 2.1(0.7) 
Total 



1.364632 5.1(0.1) 
1.364685 4.0(0.1) 
Total 



1.364807 4.7(0.5) 
1.364897 2.5(0.7) 
1.364963 4.2(0.6) 
1.365025 2.9(0.4) 
Total 



9.32(0.65) 
11.37(0.01) 
11.37(0.01) 



11.09(0.01) 
10.97(0.01) 
11.75(0.01) 
11.21(0.01) 
11.97(0.01) 



12.46(0.01) 
13.46(0.01) 
13.50(0.01) 



13.00(0.00) 
14.31(0.03) 
13.41(0.09) 
14.16(0.18) 
14.59(0.07) 



^ 13.0 
14.80(0.16) 
14.80(0.16) 



13.73(0.08) 
14.31(0.27) 
14.46(0.47) 
15.06(0.49) 
15.23(0.34) 



Velocity range B 
13.03(0.02) ^ 11.50 

13.69(0.01) 13.07(0.01) 
13.77(0.01) 13.07(0.01) 

Velocity range C 
13.18(0.01) 12.40(0.03) 
13.60(0.02) 12.62(0.10) 
14.14(0.02) 13.40(0.04) 
13.87(0.02) 13.12(0.06) 
14.43(0.01) 13.65(0.03) 



12.82(0.06) 
13.13(0.04) 
13.30(0.03) 



12.83(0.04) 
12.72(0.07) 
12.73(0.07) 
12.71(0.07) 
13.35(0.03) 



11.70 
12.92(0.01) 
12.92(0.01) 



11.98(0.05) 
11.71(0.25) 
12.63(0.01) 
12.55(0.03) 
12.97(0.02) 



sS 12.00 
12.26(0.02) 
12.26(0.02) 



^ 12.00 
^ 12.00 
12.10(0.07) 
11.97(0.11) 
12.35(0.08) 



SS18.66 
SS18.56 
^19.07 



^18.62 
SS18.34 
SS18.57 
^18.41 
^19.26 



1.365116 5.4(0.7) 
1.365217 6.8(0.7) 
1.385305 3.9(0.4) 
Total 



Velocity range D 

11.14(0.01) 13.42(0.06) 13.87(0.08) 13.62(0.01) 12.93(0.06) 12.26(0.16) 12.43(0.02) 12.34(0.02) ^18.69 

11.38(0.01) 14.16(0.80) 14.14(0.05) 14.06(0.01) 13.25(0.04) 12.63(0.08) 12.57(0.01) ^ 12.00 sS18.78 

10.78(0.02) 13.99(0.02) 13.68(0.08) 13.83(0.02) 12.69(0.07) 12.04(0.24) 11.96(0.04) 11.96(0.05) SC18.54 

11.64(0.01) 14.43(0.43) 14.41(0.04) 14.35(0.01) 13.50(0.03) 12.85(0.07) 12.86(0.01) 12.49(0.03) SC19.31 



1.365430 5.4(0.3) 
1.365570 3.9(0.2) 
1.365667 2.5(0.6) 
1.365759 3.3(0.4) 
Total 



10.70(0.02) 
10.48(0.03) 
10.43(0.04) 
10.43(0.04) 
11.12(0.01) 



13.20(0.00) 
13.69(0.01) 
12.74(0.00) 
13.11(0.01) 
13.92(0.01) 



13.42(0.08) 
13.37(0.09) 
13.37(0.10) 
13.19(0.10) 
13.95(0.05) 



Velocity range E 
13.37(0.01) 
13.52(0.01) 
13.04(0.01) 
13.10(0.01) 
13.90(0.01) 



12.51(0.01) 
12.52(0.01) 
12.19(0.03) 
12.24(0.02) 
12.99(0.01) 



12.28(0.15) 
12.11(0.19) 

^ 12.00 
11.64(0.54) 
12.56(0.14) 



11.95(0.06) 
^ 11.70 
11.70 
sS 11.70 



sS 12.00 
^ 12.00 
12.00 
sS 12.00 



SS18.69 
^18.55 
SS18.35 
SC18.46 
^19.29 



1.365906 
1.366008 
1.386197 
1.386330 



8.0(0.2) 
1.8(0.1) 
3.5(0.4) 
1.9(0.7) 



^ 9.47 
10.08(0.09) 
^ 9.47 
^ 9.47 



12.21(0.01) 
12.13(0.01) 
11.48(0.02) 
11.23(0.02) 



13.38(0.11) 
12.72(0.31) 
sS 13.0 
!£ 13.0 



Velocity range F 
12.65(0.02) 
12.30(0.05) 

^ 12.00 

< 12.00 



^ 11.50 

^ 11.50 

^ 11.50 

< 11.50 



^ 12.00 

^ 12.00 

^ 12.00 

< 12.00 



^ 11.70 
11.70 
^ 11.70 
sS 11.70 



^ 12.00 

sS 12.00 

^ 12.00 

SC 12.00 



^18.86 
18.21 
^18.49 
SC18.23 



Velocity ranges are defined in Fig. 1. 

^ Column densities are given in logarithmic units. 

^ upper limit on A^(H l) obtained from the GMRT spectrum assuming Tg — 100 K; for higher Tg this should be multiplied by (Ts(K)/100) 
b blended with other lines 



terval v and v+dv for an optically thin cloud is related to 
the optical depth, r(v), and to the spin temperature (Ts) by 



(1) 



iV(H I) = 1.835 X 10"° Ts / t{v) dv cm 



If we assume the absorption line to have a perfect Gaussian 
profile with a peak optical depth rp , then the above equation 
will become, 



iV(H I) = 1.93 X 10^* TP Ts Ay cm" 



(2) 



where AV = 2Vhi2x6 is the FWHM of the Gaussian profile 
with b the Doppler parameter. 

Note that the above equations assume the absorbing 
cloud to be homogeneous and to cover the background 
source completely. PKS 0237—233 is a very compact ra- 
dio source that is not resolved in our GMRT observations. 
VLB A observations of this source at 2.32 and 8.55 GHz have 
revealed structures at milli-arcsec scales (Fey et al. 1996). 
Most of the flux is in two components that are roughly sep- 
arated by ~ 10 milli-arcsec (see also Fomalont et al. 2000). 
These components account for all the flux at 2.32 and 8.55 
GHz implying there is no broad diffuse radio emission in this 
source. The strongest of these components has 60% of the 
flux at 2.32 GHz (Fey et al 1996) and 70% of the flux at 5 
GHz (Famalont et al 2000). It has therefore a steep spec- 
trum (/y oc and probably dominates the flux at 610 



MHz. The second brightest component has a flat spectrum 
{fi, oc i^""'^) and should correspond to the core of the radio 
source and therefore to the optical counterpart. If this is the 
case, then the optical and 610 MHz radio lines of sight are 
separated by 10 milli-arcsec. However, the spectral indices 
derived above are highly uncertain due to different beam 
sizes and component decompositions of the two above ob- 
servations. Because of this, it is unclear which of these com- 
ponents is indeed associated with the optical point source 
and the angular separation of 10 milli-arcsec should be con- 
sidered as an upper limit. One milli-arcsec a.t z = 1.35 cor- 
responds to a linear size of 8.5 pc (for a flat universe with 
Qa = 0.73, = 0.27, Ho = 71 km/s/Mpc). This means 
the optical and radio sightlines are either coincident or sep- 
arated by ~85 pc. 

Assuming that the velocity width of the H I feature is 
the same as that of metals, and using rp(3a) = 0.003 and T 
= 100 K, as expected in the cold neutral medium, we derive 
upper limits on A''(H l) in individual components. These val- 
ues are listed in last column of Table. [1] These limits should 
be scaled up by Ts(is')/100. Also for temperatures typical of 
the warm neutral medium (WNM) or warm ionized medium 
(WIM) of our Galaxy, higher b values (say 9 km s~^) should 
be considered. 

It is clear that the uncertainties on the temperature and 
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Figure 4. Top panel shows the GMRT spectrum on a veloc- 
ity scale centered at Zabs = 1-3647. The horizontal dashed lines 
correspond to the mean and la range in the normalised flux. 
The corresponding C 1AI6O6 absorption is plotted in the bottom 
panel. Vertical tick marks show the location of diflerent velocity 
components derived from the metal line profiles. 
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Figure 5. Input spectral energy distributions (SEDs) used in 
our model calculations. Continuous and dashed lines are SED 
taken from Haardt & Madau (2005) for the UV extra-galactic 
background contributed respectively by QSOs alone and QSOs 
and Lyman break galaxies. Stars indicate the Galactic UV back- 
ground SED computed by Black (1987). The vertical dashed line 
shows the Lyman Limit; the ionization potentials of and Mg" 
are shown with arrows. 



Table 2. Total column density: 



Species 


log iV(cm"^) 


N(H i) 


19.30±0.30 


N(0 I) 


>15.46 


N(Mg I) 


12.39±0.00 


N(Mg II) 


^14.98 


N(C II*) 


13.73±0.02 


N(Fe II) 


14.13±0.01 


N(Si II) 


14.89±0.01 


N(A1 II) 


^ 13.46 


N(A1 III) 


13.53±0.01 


N(Mn II) 


11.81±0.03 


N(Cr II) 


12.32±0.13 



the covering factor of the cloud in front of the radio-source 
prevents any rapid conclusion. To go further in the analysis 
of this system, in the following we construct photoioniza- 
tion models and constrain them using information from the 
21 cm, Lya , C I, C I* and C 11* absorption lines. 



4 MODEL CALCULATIONS 

As the total Af(H l) in this system is less than 10'^° cm~^, it 
is most likely that ionization correction is important when 
deriving metallicity for most of the species. As an exception 
to this, O I and H I are coupled through charge-exchange 
reaction and for iV(H i)^ 10^'' cm~^ iV(0 i)/N{}i i) di- 
rectly reflects [O/H] (see Viegas, 1995 and results presented 
below). As O 1AI3O2 is possibly saturated in the central 
components (see Fig. [1]) we can derive only a lower limit on 
N{0 l). We find [0/H]> -0.33 if we consider log 7V(H l)~ 



19.3. High metallicity is also inferred from Mg 11 and Si 11 
(see Table. [2]). Clearly the gas under consideration must have 
gone through (or is going through) a prolonged period of 
vigorous star-formation activity. 

To correct for ionization we use model calculations per- 
formed with the photoionization-simulation code CLOUDY 
(Ferland et al. 1998). Keeping in mind the preferred value 
log A'^(H i) — 19.3±0.30, we consider three total H i column 
densities, log N{R l) = 19.0, 19.5 and 20.0, at which the 
calculation is stopped. Relative metal abundances are con- 
sidered solar and the total metallicity is Z = 0.5 Zq. This is 
motivated by the above inferred oxygen metallicity. 

Most of the model calculations discussed below use a 
ionizing radiation spectrum that is a combination of the 
meta-galactic UV background at z = 1.35 contributed by 
QSOs and galaxies as computed by Haardt & Madau (2005) 
(dashed SED in Fig. [5]) and the Cosmic Microwave Back- 
ground with TcMBR = 6.4 K. Relative to quasars, galaxies 
contribute approximately an order of magnitude higher flux 
in the energy ranges corresponding to the ionization poten- 
tials of C'' and Mg'' (See also Wolfe et al. 2004). We also 
consider models with additional local radiation (continuous 
curve with points in Fig. [5]). For this purpose we use the 
Black (1987) galactic background radiation after removing 
the H I ionizing photons as usually done in the case of ISM 
simulations (see e.g. Draine & Bertoldi, 1996; Shaw et al. 
2006). We consider dust to metal ratio 1/10 of that observed 
in the Galactic ISM. 

First we mainly focus our attention on models with 
meta-galactic radiation field contributed by QSOs and 
galaxies. For a given values of total A'^(H l) and metallic- 
ity (Z) we generate grids of models varying the density in 
the range, —3 ^ log nH(cm~^) ^ 1.5. The results of our 
model calculations are summarized in Figs. |6] and [T] From 
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Figure 6. Results of photoionization models. The solid, short-dashed and long-dashed curves are for, respectively, log A'^(H l) = 20.0, 
19.5 and 19.0. Metallicity is 0.5 Zq. Constraints from the observed ratios of column densities integrated along the line of sight are shown 
in each panel as boxes or arrows. The arrows in panel (b) indicate the upper limits on njj obtained from the upper limit on the 21-cm 
optical depth assuming the spin temperature to be equal to the kinetic temperature and complete coverage of the background radio 
source. 



these figures it can be seen tiiat, for the range of A'^(H l) 
considered here, most of the constraints are consistent with 
0^ log[nH (cm^"')] ^ 1. This density range corresponds to 
the gas being part of a CNM phase (T < 100 K) only for 
the model with log N{R l) = 20 (see panel (a) of Fig. IHl). In 
the other two models, that have A'^(H l) consistent with the 
observed range, the gas has the characteristics of a WNM 
phase or is in the thermally unstable region. Therefore, it 
seems most likely that the C I absorption is originating either 
from WNM or WIM. The N{Ug l) /iV(Mg ll) ratio predicted 
by these models is larger than the observed ratio. This can 
be explained by the large uncertainties in N{Mg ll) due to 
heavy saturation of the absorption lines (see Fig. 1). 

The models allow us to estimate metallicity and deple- 
tion. From Fig.[7]it is clear that for —1 ^ log nH(cm~^) ^ 1 
the ionization corrections for H I, O I, Mg ll and Si ll 
are negligible. Thus these ionization states can be used di- 



rectly to infer the overall metallicity of the system. For 
N{U l)~ 10^^'^ cm~^ the observed column density of Si ll 
and the lower limits on O I and Mg ll column densities are 
consistent with Z = 0.56 Zq. For such a metallicity the pre- 
dicted column densities of Fe ll and Mn ll are about an order 
of magnitude higher than the observed column densities. 
The lower gas-phase abundance of Iron co-production ele- 
ments with respect to the a— process elements is confirmed 
by the measured low column density of Cr ll and the absence 
of Ni II absorption. If this difference is due to depletion onto 
dust, then predictions for the Zn ll column density should 
match the observed value. This is not the case and Zn ll is 
overpredicted by the models. This argues toward an excess 
of a— process elements. 

Before drawing strong conclusions on the physical state 
of the gas, size of the absorbing region and the metal en- 
richment scenario it may be interesting to investigate the 
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Figure 7. Results of photoionization models. Various ion column densities predicted by the models are plotted as a function of hydrogen 
density Symbols are as in Fig. E] /(H l) = Af(H l)/N(H). 



physical state of C I sub-components. While single cloud 
models give the average property of the gas, modelling indi- 
vidual velocity ranges will allow us to investigate variations 
of physical conditions across the profile. In the following 
Section we therefore discuss properties of the different C I 
velocity ranges we have identified before (see Fig. [1} . In the 
corresponding models we have varied A''(H l) from 10^* to 
10^^ cm-^ 



4.1 Discussion on the different velocity ranges 
4-1.1 Velocity range A 

This velocity range has three well resolved components and 
an additional weak component in the red wing. C I ab- 
sorption is possibly present in the three main components 
with little C I* absorption. We find that the upper limit on 



N{C I*) is consistent with log nH(cm~^) ^ 0.0 for mod- 
els with log A''(II l)^19. C II* absorption is blended with a 
strong Ly-a line at 2abs= 1-5. Therefore, there is no strong 
constraint on the density for this velocity range. This ve- 
locity range has log N{Mg l)/7V(Mg ll) ~ —2.4 and log 
N{0 i)/N{Si ll) = 0.06±0.05. In our model calculations 
these ratios are consistent with log N{H l) ^18.5 and 
—3 ^logriH (cm~'^) ^ —2.0. Such low densities are needed to 
explain that iV(0 l)/iV(Si ll) is different from [0/Si]. We also 
find that the observed log N{Al lll)/7V(Al ll) = -0.22±0.11 
ratio is consistent with low 7V(H l) and low density. The 
presence of possible additional local radiation below the Ly- 
man limit from the local stellar population would reproduce 
the observed A'^(Mg l)/A'^(Mg ll) ratio at higher densities 
(Fig. . In any case models for this gas require large ion- 
ization corrections for H I (see Fig. [7}. We therefore do not 
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Figure 8. Profiles of high-ionization lines detected at ^abs = 
1.3647. Tick marks indicate the location of the different velocity 
components seen in the low-ionization profiles. Vertical dashed- 
lines mark the location of components that are detected only in 
the high ionization species. 



Table 3. Model results for the velocity range B 



Species 


log N 


log 


Ar(H i) 


= 19.0- 


18.5 




(Observed) 


Z 


= 0.1 Zq 


Z = 0.5 Z 


C I 


12.92±0.01 


11 


99 


-12.14 


12.66 


-12.84 


C I* 


12.56±0.02 


11 


28 


-11.60 


11.95 


-12.30 


C II* 


13.30±0.03 


12 


14 


-12.38 


12.81 


-13.06 


O i" 


14.23±0.05 


14 


24 


-14.72 


14.94 


-15.42 


Mgi 


11.37±0.01 


11 


50 


-11.74 


12.08 


-12.31 


Mg II" 


13.50±0.01 


13 


40 


-13.77 


14.06 


-14.44 


Al II 


12.22±0.01 


12 


41 


-12.85 


13.05 


-13.47 


Si II 


13.77±0.01 


13 


48 


-13.93 


14.30 


-14.45 


Fe II 


13.07±0.01 


13 


30 


-13.70 


13.99 


-14.42 


Mn II 


10.93±0.07 


11 


25 


-11.67 


11.98 


-12.39 


Zn II 


10.61±0.20 


10 


34 


-10.66 


10.99 


-11.30 


Cr II 


^ 11.30'' 


11 


37 


-11.80 


12.16 


-12.56 



" probably saturated; 3(7 upper limit. 

attempt to further infer metallicities and depletion pattern 
in this velocity range. 



4-1-2 Velocity range B 

This velocity range has two components (one main compo- 
nent and a weak blend in the blue wing) and is associated 
with the strongest C I and C ll* absorptions. From Fig. [l] 
it appears that the O l absorption is not heavily saturated. 
This velocity range contributes about ~20% of the total 
A''(0 l). This means that about 20% of the total iV(H l) 
is in this velocity range and A'^(H l) is probably not much 
higher than lO" cm"^ The observed log N{0 i)/A^(Si ll) = 
1.03±0.13 and log iV(0 i)/A^(Mg ll) = 1.30±0.16 ratios are 
consistent with solar metallicity ratios [0/Si]© = 1.13 and 



[O/Mg]0 = 1.11. As all these species are a— process elements 
it is most likely that ionization corrections are negligible. 

The measured A'^(C l*)/A'^(C l) is consistent with 
-1.0 ^ log[nH(cm~^)] ^ -0.6 for 18 log [iV(H l) 
(cm~'^)] ^ 19.5. The temperature of the gas is higher than 
1000 K for the range in metallicity and log A'^(H l) considered 
here. The gas pressure in this component is consistent with 
3.0 ^ log[P/k(cm-^K)] < 3.6. The gas is therefore warm 
with a neutral hydrogen fraction ^0.7. 

The observed A''(C ll*)/A''(Si ll) ratio is much higher 
than that predicted by the models in the range of nn de- 
rived from N{C l*)/N{C l). The difference is higher for 
higher A^(H l) models. In Table. [3] we give the column den- 
sities of different species predicted in our model for — 1.0 ^ 
log[nH(cm*^)] < -0.6, 18.5 log [iV(H l) (cm^^^)] < 19 
and two different assumed values of metallicities (Z = 0.1 
and 0.5 Zq). From this table it is clear that the observed 
column densities of Si ll, Mg ll, Zn ll and O I are well repro- 
duced by the models with Z = 0.1 Zq. However, the observed 
column densities of C I, C I* and C ll* are higher than in 
the models. Also the well measured column densities of Fe ll 
and Mn ll are a factor of two below the model predictions. 
Models with Z ^ 0.5 Zq are needed to explain A'^(C l) and 
N{C II*). These models will require a depletion of O, Mg, 
Al and Si by up to 0.6 dex and Fe, Mn and Cr by up to 1.0 
dex. 

The observed ratio, log 7V(A1 lll)/iV(Al ll) = -0.15 ± 
0.17, is similar to what is observed for component A. This 
clearly suggests the presence of low-density ionized gas (with 
low column density) over the same velocity range as the 
C I gas. This is confirmed by the presence of strong Si iv 
and C IV absorption lines with similar velocity structure 
(see Fig. |5}. It is interesting to note that the observed 
A'^(Mg l)/A''(Mg ll) ratio is explained by the models dis- 
cussed here. This means that there is no excess of UV radia- 
tion with energies less than 13.6 eV in these velocity compo- 
nents (see Fig. [9} . In conclusion this velocity range seems to 
be well apart from any local star-forming region, has metal 
enrichment between 0.1 and 0.5 Zq and reasonable depletion 
of iron co-production elements. 



4-1-3 Velocity range C 

This velocity range has ^ 60% of the total O I column den- 
sity distributed in 4 distinct velocity components. Only two 
of them show detectable C I* absorption and other two com- 
ponents are mainly in the wings of the absorption. Clearly 
most of the neutral hydrogen in the system must originate 
from this velocity range. The total column densities of C I, 
C I* and C ll* in this velocity range are consistent with that 
seen for velocity range B. This means that the average pres- 
sure in both the velocity ranges are more or less consistent 
with one another. The measured A''(C l'')/N{C l) ratio is 
consistent with —1.0 ^ log[nH(cm~^)] —0.2 for the range 
of N{}i l) considered in our models. We also find that the 
models with 18.5 ^log \N{Y{ l)(cm"^)] ^ 19 wiU also re- 
produce the observed A'^(C ll*)/A''(Si ll) ratio. This means 
that we do not require enhancement of C with respect to Si 
and Mg as in component B. The gas pressure in this veloc- 
ity range is 3.0 ^ log[P/k(cm~'^K)] ^ 3.8. Like in velocity 
range B, the C I absorption from this component comes from 
a WIM or WNM. 
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Figure 9. Results of photoionization models with an excess of local radiation in addition to the meta-galactic UV background radiation. 
Solid and dashed lines are for, respectively, log N(ll l) = 19.0 and 18.5. Thick and thin curves are for models with and without local 
radiation added to the meta-galactic background radiation. 



O I and Mg ll absorptions are saturated and the derived 
column densities of Si ll, Al ll, Fe ll, Mn ll and Cr ll are 
higher than in velocity range B. The column densities of 
C I, C 11*, O I, Si 11 and Mg ll are approximately consistent 
with Z 0.5 Zq (see Table|4]|. However the observed column 
densities of Fe ll, Mn ll and Cr ll and the 3 a upper limits 
on Ca 11 and Ni ll are consistent with the iron co-production 
elements being under abundant by ~ 0.4 dex. This depletion 
of iron co-production elements with respect to a— elements 
is consistent with what is seen in the velocity range B. 

As we have only a lower limit on the Mg ll column 
density, the actual A'^(Mg l)/A'^(Mg ll) ratio may be smaller 
than that seen in velocity range B. From panel (e) in Fig. [9] 
it is clear that for log A'^(H l) = 19, the observed log 
A'^(Mg l)/A'^(Mg ll) ratio is consistent with a UV radiation 
field as high as the galactic mean UV field below 13.6 eV. 



However, the observed A'^(Mg l) and A''(C l) are inconsistent 
with much larger radiation fields (of the order of 10 times the 
Galactic field) as required for the H2 components in high-z 
DLAs. 

Like in the case of velocity range B we notice that 
the Al III column density is similar to that of Al 11 (log 
N{Al lll)/A''(Al 11) = -0.04 ± 0.02). Also very strong Si iv 
and C IV absorptions are present with similar velocity struc- 
ture as the low ionization lines (see Fig. [S| . This again sug- 
gests that highly ionized gas is co-spatial with the C I com- 
ponent in this velocity range as well. 



4-1-4 Velocity range D 

The O I absorption in this velocity range is not saturated. 
However, Mg 11 lines are saturated (see Fig.[T}. Only one out 
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Table 4. Model Results for the velocity range C 



Species N log 7V(H i) = 19.0-18.5 





(Observed) 


Z 


= 0.1 Zq 


Z : 


= 0.5 Z 


' 


C I 


12.97±0.02 


12, 


,04- 


-12, 


,14 


12, 


.66- 


12, 


.88 


C I* 


12.65±0.10 


11, 


,28- 


11, 


,98 


11, 


.94- 


12, 


.54 


C II* 


13.35±0.03 


12, 


.14- 


-12, 


,47 


12, 


.81- 


-13, 


.13 


O I" 


^ 15.23 


14, 


.22- 


-14, 


,72 


14, 


.92- 


15, 


.42 


Mg I 


11.97±0.01 


11, 


.43- 


-11, 


,74 


11, 


.99- 


-12, 


.31 


Mg II" 


^ 14.59 


13, 


,30- 


-13, 


,77 


13, 


.98- 


-14, 


.44 


Al II 


13.08±0.01 


12, 


.29- 


-12, 


.85 


12, 


.95- 


-13, 


.47 


Si II 


14.43±0.01 


13, 


.37- 


-13, 


.93 


14, 


.04- 


-14, 


.45 


Ca II 


< 12.00'' 










12, 


.45- 


12, 


.68 


Mn II 


11.53±0.04 


11, 


.20- 


11, 


.67 


11, 


.91- 


12, 


.39 


Fe II 


13.65±0.03 


13, 


.23- 


-13, 


.70 


13, 


.90- 


14, 


.42 


Ni II 


^ 12.95'' 










12, 


.76- 


13, 


.34 


Zn II 


< 11.10'' 


10, 


.28- 


10, 


.66 


10, 


.94- 


11, 


.30 


Cr II 


11.63±0.14 


11, 


.34- 


-11, 


.80 


12, 


,11- 


12, 


.56 



" probably saturated; '' 3a upper limit. 



Table 5. Model results for range D 



Species 


N 


log N{U i) = 18.0 




(Observed) 


Z = 0.5 Z Q 


C I 


12.86±0.01 


12.47-12.53 


C I* 


12.34±0.02 


11.65-12.52 


C II* 


12.85±0.07 


12.60-12.92 


O I" 


14.41±0.04 


14.43 


Mg I 


11.64±0.01 


11.75-12.17 


Mg II" 


^ 14.43 


13.70-14.25 


Al II 


13.08±0.01 


12.67-13.38 


Si II 


14.35±0.01 


13.76-14.45 


Fe II 


13.50±0.03 


13.59-13.97 



° probably saturated; '' 3a upper limit. 



of the three components show detectable C I* absorption. All 
components have N{0 l)~ A^(Si ll) within 0.2 dex. For mod- 
els considered here, N{C l*)/N{C l) and iV(C ll')/iV(Si ll) 
are consistent with, —1.2 log[nH(cm~"^)] ^ 0.0 and 
2.7 < log[P/k(cm-3K)] < 3.8. 

Models predict log iV(0 i)/iV(Si ll)^ 0.9 for log iV(H l) 
> 18.5. Therefore observations are either consistent with 
log A'^(H l)<18.5 or require under-abundance of O with re- 
spect to Si. It can be seen from Table [5] that most of the 
column densities of singly ionized species can be reproduced 
by models with log iV(H l) = 18 and Z = 0.5 Zq. 

4.2 Discussion 

Rao & Turnshek (2000) have shown that DLAs can be pre- 
selected based on the equivalent widths of Mg ll, Fe ll and 
Mg I absorption lines. Specifically, they have found that 50% 
of the absorbers with rest equivalent width of Fe llA2600 
and Mg llA2796 greater than 0.5 A were confirmed DLAs. 
If they impose the condition iy(Mg 0.5 A then all the 
systems are confirmed DLAs. It is also clear from Fig. 3.4 
of Lane's PhD thesis (2000) that 50% of these systems will 
be detected by 21-cm absorption. The inferred low A'^(H l) 
from lUE data (log A''(H i) ~ 19.3) and the absence of 21-cm 
absorption at Zahs= 1.3647 towards PKS 0237-233 are both 



consistent with the above findings as VKr(Mg l)~0.3A in this 
system which is clearly a sub-DLA. It is therefore somewhat 
surprising to have found a wide-spread C I absorption as- 
sociated with the system suggesting high metallicity. Wide- 
spread (i.e AV ^ 100 km s~^) C i absorption is seen at z^hB= 
1.973 towards Q 0013—004 (Petitjean, Srianand & Ledoux 
2002) and 2abs= 1.962 towards Q 0551-366 (Ledoux, Sri- 
anand & Petitjean 2002). Molecular hydrogen is detected 
in these systems although several of the C I components do 
not have associated H2 absorption. C i absorption associated 
with II2 have been shown to originate mainly from CNM gas 
(Srianand et al. 2005). However the nature of the gas in the 
components that do not have associated II2 absorption is 
still unclear. 

In the case of the Zabs^ 1.365 system toward 
PKS 0237—233 a simple one cloud model can explain the to- 
tal column densities. Under this assumption of single cloud, 
it can be shown that the gas producing C l absorption is 
most unlikely to be a CNM gas if the actual A'^(H l) is close 
to that we measure from the lUE spectrum. Detailed mod- 
elling of individual velocity ranges are also consistent with 
the gas responsible for the C I absorption in this system 
being either WNM or WIM. In addition we also infer ion- 
ization and chemical differences between individual velocity 
ranges. 

Our GMRT observations give a good constraint on 
A''(II l) /Ts if we assume complete coverage of the background 
radio-source by the absorbing gas. However, as pointed out 
before, the background radio source shows structures at 
milli-arcsec scale with a separation between the two bright- 
est components of 85 pc at the absorption redshift. It is un- 
clear which of these components coincides with the optical 
point source. If the brightest component that has most of the 
flux in the redshifted 21-cm range coincides with the optical 
source then we would have detected 21-cm absorption if the 
C I gas originates from the CNM. The non-detection would 
be a definite argument for the gas being warm. However, 
if the weaker radio component coincides with the optical 
source, then even if the C I absorption originates from CNM 
gas our GMRT observations would not be able to detect the 
corresponding 21-cm absorption as the expected size of the 
cloud is smaller than the radio source separation. It is known 
that strong Mg 11 systems (with Wr >0.3 A) show coherent 
absorption over more than ~ 2 kpc (Petitjean et al. 2000, 
Ellison et al. 2004). The fact that we find 21 distinct com- 
ponents along the line of sight to the optical point source 
means that it is likely that the overall system covers both 
the radio components. In that case our GMRT observations 
are consistent with no CNM with log A''(H l)^ 19 along the 
radio source as well. Thus, it seems most likely that the re- 
gion probed by the Mg 11 absorption system at Zabs= 1.365 
is composed mostly of a warm neutral medium and/or warm 
ionized medium. 

The inferred metallicities suggest that the gas has gone 
through an active star-formation period. We do find the sig- 
nature of strong Q— process element enhancement with re- 
spect to Fe co-production elements. However, as we have 
not detected the Zn 11 absorption it is not clear, although 
unlikely, whether the difference could be attributed to dust- 
depletion. Photo-ionization models for the individual com- 
ponents are consistent with ionization by meta-galactic UV 
background. However, inconsistent with very large local ra- 
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diation field as one sees in H2 components in high-z DLAs. 
This suggests that the system is, at the time of observation, 
in a low star-formation state. 

Out of the 50 Mg 11 systems studied in the ESO large 
programme 'The Cosmic Evolution of the IGM', only two 
systems show C I absorption (i.e 2abs= 1.365 and 1.6724 
towards PKS 0237-233). The only other sub-DLA towards 
radio quiet QSO HE 0001-2340 does not show C I and other 
fine-structure lines (See Richter et al. 2005). The inferred 
N{}i l), metallicity, and depletion pattern at 2abs= 1.365 
toward PKS 0237—233 are similar to what is observed at 
0abs= 2.139 toward Tol 1037-2703. The latter system also 
shows detectable C l* and C 11* absorption lines spread over 
~ 70 km s~^and distributed in two distinct velocity com- 
ponents without detectable H2 absorption. Interestingly the 
only other system that shows C I absorption in the large 
programme data is also found towards PKS 0237—2703. 
This line of sight is famous for the presence of a super- 
cluster of C IV absorption lines. Thus, our study provides 
yet another motivation to go for further deep imaging and 
follow-up spectroscopic observations of the field (see Aragon- 
Salamanca et al. 1994). 
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